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Anomalous Hall Effect in Ferromagnetic Metals: 
Role of Phonons at Finite Temperature 
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The anomalous Hall effect in a multiband tight-binding model is numerically studied taking into account both elastic 
scattering by disorder and inelastic scattering by the electron-phonon interaction. The Hall conductivity is obtained as 
a function of temperature T, inelastic scattering rate 7, chemical potential fi, and impurity concentration Xi^p. We 



find that the new scaling law holds over a wide range of these parameters; —ax 
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\Gxx + b, with 



Gliv (cr^^o) being the conductivity tensor (with only elastic scattering), which corresponds to the recent experimental 
observation [Phys. Rev. Lett. 103 (2009) 087206]. The condition of this scaling is examined. Also, it is found that the 
intrinsic mechanism depends on temperature under a resonance condition. 
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The anomalous Hall effect (AHE) is a phenomenon where 
the transverse conductivity arises owing to the relativistic 
spin-orbit interaction in the presence of the spontaneous mag- 
netization.' This effect has attracted much attention because 
it is closely relevant to many fundamental questions on quan- 
tum transport in solids. Historically, the first theory of AHE 
was proposed by Karplus and Luttinger (KL) who considered 
the matrix elements of a current operator between different 
bands and discovered the so-called anomalous velocity.^ In 
their original paper, they did not mention the disorder effect 
and were criticized by Smit who stressed that the accelerat- 
ing electric field is balanced by that of the impurity potential 
in the steady state. He instead proposed the skew scatter- 
ing mechanism, where the transition probability for the scat- 
tering fc — >^ fc' is different from that for k' — s> k. Namely, 
the detailed balance is violated in a scattering process, lead- 
ing to the transverse current to the electric field. Later, Berger 
proposed another mechanism called the side jump, where a 
transverse shift in the electron trajectory occurs during the 
scattering.^' ^ These two are collectively dubbed the extrinsic 
mechanisms, while that of KL the intrinsic mechanism. For 
a long time, it has been considered that the extrinsic mech- 
anisms are dominant, while the intrinsic mechanism is neg- 
ligible. ' Experimentally, the Hall resistivity is analyzed 
using the expression 



tot 
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RqH + Pxy 



(1) 



where i?o is the ordinary Hall coefficient, while the second 
term p^y is the anomalous contribution due to the spontaneous 
magnetization Ms and is often empirically written as AttRsMs 
with the anomalous Hall coefficient R^. Furthermore, the re- 
sistivity (pxx) dependence of p^y has been analyzed using 



Pxy = apxx + f3p: 



xx ' 



(2) 



where the first term corresponds to the skew scattering, while 
the second one coresponds to the side jump. However, note 
that the intrinsic contribution could also contribute to the sec- 
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ond term. 

Recent development has been achieved by theoretical re- 
formulation based on the Berry curvature in the momentum 
space and by first-principles band calculations,'^"' ' which are 
combined to reveal the topological nature of the AHE.' The 
disorder effect on the AHE has been theoretically studied to 
reveal the scaling behavior,'^' '^ which is consistent with the 
systematic experimental results.'"'^ However, most theoret- 
ical analyses have been limited to the zero temperature limit 
where inelastic scattering is neglected. The AHE at finite tem- 
perature has already been studied,'^ '^ but remains controver- 
sial. Therefore, a more systematic study including inelastic 
scattering at finite temperature is highly desired. 

Experimentally, it has been found that the skew scattering 
contribution is rapidly suppressed as temperature increases, 
and the Hall conductivity axy reaches a rather steady value. ' 
More explicitly, an empirical formula for cTxy has been pro- 
posed by Tian et al. '^ as 



^xy 



{T,d)^plfid)aUT,d) + b, 



(3) 



plfid) = a(d)a-i(T ^ 0,d)+P{d)a^^{T ^ 0,d), (4) 

where crxx{T — > 0, d) is the residual conductivity at the low- 
est temperature experimentally accessible. They showed that 
this expression fits well the data in Fe thin films with vary- 
ing temperature T and the thickness d. The remarkable ob- 
servation in ref. 18 is that b is independent of both T and d, 
which is identified as —cr^xy- Note, however, that the temper- 
ature dependence of h was reported for Ni thin films. '^ The 
ideas underlying eq. (3) are that the roles of elastic and in- 
elastic scatterings are different, and that only the former con- 
tributes to the extrinsic mechanisms, while the latter rapidly 
suppresses them. On the other hand, the intrinsic contribution 
a^xy is expected to be much more robust. Although eqs. (3) 
and (4) seem to work well experimentally, their justification 
and theoretical explanation are still lacking. 

In this paper, we study both the intrinsic and extrinsic 
mechanisms in a simple model taking into account both elas- 
tic and inelastic scatterings at finite temperature T much 
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Struct a fully spin-polarized multiband tight-binding model on 
the square lattice: 
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Conduction electrons are described by the creation and anni- 
hilation operators c\ (c; ) . Each impurity has the s and — ip^ 
orbitals, which are described by s|(si) and pj(pi), respec- 
tively. The phase factor e^'^'^ originates from the spin-orbit- 
coupled p^ — ipy orbital, where Oij is the angle from the impu- 
rity at the jth site to the electron at the ith site measured from 
the X axis. For simplicity, the transfer integrals between im- 
purities are neglected. Also, the impurity concentration Ximp 
is changed from zero to unity. Therefore, this model connects 
two different limits. One is the periodic three-band model in 
the dense limit Ximp — 1, which exhibits the intrinsic mecha- 
nism due to the existence of the p^ — ip^ orbital. ' The other is 
the single-impurity model in the dilute limit x-^^p <^ 1, which 
exhibits the skew scattering due to the interference between 
the s and p orbitals. Generally, two orbitals with the azimuthal 
quantum numbers / and I + 1 are necessary for the skew scat- 
tering. ^""^^ Below, we choose = 1, — ei =£2 = 1, and 
Vi =¥2 = 0.1, which leads to a reasonable Hall angle on the 
order of 10~^. 

Therefore, our model is the minimal one containing all the 
essential features of the AHE at finite temperature. We will 
consider only the electron-phonon interaction as the source of 
inelastic scattering, and neglect the electron-magnon interac- 
tion. This is justified when T, 0d <C Tq with the Debye tem- 
perature 6d. This condition is satisfied for Fe (where the new 
scaling laws eqs. (3) and (4) have been proposed'^) and Co 
with Tc > lOOOK, but is marginal for Ni with Tq = 627K, 
and is not justified for materials with lower Tq. 

The conductivity tensor is numerically calculated by the 
Kubo formula, 

27r ^ /(U) - /(Cn) {m\J^^\n){n\r\m) 
C't^u = > ^ j — ' (") 



mn 



for each impurity configuration. A set of conductivities are 
averaged over 9600 configurations for < a;iinp < 0.3, 6400 
for 0.4 < Xin,p < 0.6, and 3200 for 0.7 < Xi^.p < 1. 
Here /(^) = (e^^-^ + 1)^^ is the Fermi distribution func- 
tion with the energy ^ measured from the chemical potential 
/i and finite temperature T. The inelastic scattering rate 7 is 
phenomenologically introduced as the imaginary part of the 
self-energy, which in reality depends on temperature and fre- 
quency. The effect of the frequency dependence of the self- 
energy will be discussed later. The factor 2tt is due to the 
conductivity unit of e'^/h = (27r)~\ and iV^ = 30 x 30 
is the number of sites. The eigenstate \n) is represented on 
the real-space basis because we consider the general impurity 
concentration Xi^^p- Our numerical calculations can fully take 
into account the impurity effects for a given impurity con- 
figuration because we employ the numerical diagonalization 
and the Kubo formula itself is exact. The ab initio Korringa- 
Kohn-Rostoker method combined with the coherent potential 
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Fig. 1. (Color online) Calculated Hall conductivity —crxy (a) as a func- 
tion of inelastic scattering strength 7 for T = 0.1 and (b) as a function of 
temperature T for 7 = 0.1. Xjmp = 0.1, 0.5, 1.0 and fi = 62 = 1, 1.2. 



approximation was employed to investigate the AHE,^^ but it 
is applicable only to elastic scattering and zero temperature. 
On the other hand, our method can reveal a unified picture of 
the AHE at finite temperature with the inelastic scattering rate 
7- 

We can obtain the Hall conductivity as a function of T, 7, 
and Ximp. In Fig. 1(a), we can find the strong 7 dependence 
of —(Txy in the dilute regime, and the weak 7 dependence in 
the dense regime. The temperature dependence at 11 = €2 is 
strong, but is weak in other cases, as seen in Fig. 1(b). Thus, 
the calculated Hall conductivity shows different dependences 
on many parameters, and is analyzed below. First, we can find 
the 7 dependence of axx in Figs. 2(a) and 2(b), which is fitted 
by (t-HT, 7) = f^^io(r)(7 + 7o)/7o, i.e., the Matthiessen's 
rule holds. Thus, 70 obtained by fitting can be assigned to the 
elastic scattering rate due to disorder, while 7 is assigned to 
the inelastic scattering rate. As shown in Fig. 2(c), —(Txy is a 
linear function of cr^^ by changing 7 in the dilute, intermedi- 
ate, and dense regions. Namely, the relation 



^..(r,7) = Plfo(T)cTlJT,^) + b{T) 



(7) 



holds for any impurity concentration < Ximp < 1. Here, the 
subscript does not mean the T -> limit, but the 7 — 
limit. The temperature dependences of p'^q{T) and b{T) 
come from the Fermi distribution function only. In real ex- 
periments, the temperature dependence of —dxy (T, 7) mainly 
comes from that of 7 = j{T). Therefore, eq. (7) almost corre- 
sponds to eq. (3) experimentally proposed in ref. 18. By com- 
bining eq. (7) and the Matthiessen's rule, we can conclude 
that the intrinsic contribution ~<Txy — b{T) is robust against 
inelastic scattering, while the extrinsic contribution is rapidly 
suppressed as -a^f = plfo'^lxo/ (7/70 + 1)^. The extrinsic 
contribution in the first term can be further separated into the 
skew scattering and side jump contributions using 

P7y,m = a{T)^^UT) + P{T)c7-UT), (8) 
corresponding to the second scaling law eq. (4). As shown 
in Fig. 2(d), P^'^qCtIxO is ^ linear function of axxo in the di- 
lute and dense regions, respectively, but not in the interme- 
diate region because higher-order perturbations with respect 
to disorder potential are relevant and the band structure is 
ill-defined. In Fig. 3, the separated contributions of the in- 
trinsic mechanism, skew scattering, and side jump are plot- 
ted as functions of the chemical potential. In Fig. 3(a) for 
a^imn = 0.1 in the dilute limit, the skew scattering contribu- 
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a{T)al^{T, j) / axxoiT) is almost dominant 
ei and €2- This is consistent with the fact that the 
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Fig. 4. (Color online) (a) Frequency dependence of the imaginaiy pail of 
the self-energy, V{yj) = ~3S(aj). (b) Chemical potential dependences 
of the Hall conductivities calculated using eqs. (9) (points) and (6) (lines), 
xjnip = 0.1 and ©d = 0.1. 



Fig. 2. (Color online) (a) Longitudinal resistivity pxx as a function of 7 
and (b) its magnified image for Ximp = 0.1,0.3,0.5,0.7,0.9. (c) Hall 
conductivity —dxy as a function of cr^^ by changing 7, which is based on 
eq. (7). (d) Plot for p'^^yO^xxO "^^0 obtained using eq. (8). T = 0.3 and 
/i = £2 = 1. 
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Fig. 3. (Color online) Intrinsic (red filled circle), skew scattering (green 
open triangle), and side jump (blue star) contributions extracted from the total 
Hall conductivity (black filled square) as functions of the chemical potential 
fi.T = 0.1, 7 = 0.05, and (a) Xi^p = 0.1 and (b) Ximp = 0.9. 



orbitals. On the other hand, in Fig. 3(b) for ajjmp = 0.9 in the 
dense limit, the intrinsic mechanism is dominant. Note that 
the s orbital of an impurity does not contribute to the Hall 
conductivity in this limit. 

Now, we examine the validity of eq. (6) with the constant 7 
in the above calculation. The imaginary part of the self-energy 
r(cij) = — depends on frequency, and the vertex cor- 
rection exists. According to the Migdal's theorem,^^ the ver- 
tex correction is on the order of (9(6d / Ep) ~ 10^^ and neg- 
ligible. As for the self-energy, we have calculated the lowest 
order as shown in Fig. 4(a). It is found that the frequency de- 
pendence of r(a;) is appreciable there only when oj, T < 9d; 
it is negligible otherwise. This frequency dependence is taken 
into account using 

X / dcoAM I duo'A,,[u')p^^l^J^, (9) 

J J (uj-iv' + idy 

where ^m(w) = tt^'^T{lo) / {{uj — + ^^(a;)) is the spec- 
tral function and 5 is infinitesimal. 

Since the effect of inelastic scattering is mainly on the 
extrinsic mechanisms, we choose the case of Ximp = 0.1 
where the skew scattering has a large contribution. Figure 4(b) 



(6). The former is found to be approximated by the latter 
with 7/2 = T{uj — 0) rather than T{lu 00) down to 
T = 0.66d = 0.06, which is the lowest temperature at 
which we could obtain sufficiently accurate data of numerical 
calculation. This suggests that the Fermi surface term, rather 
than the Fermi sea term, is dominant in the extrinsic mech- 
anisms similarly to the longitudinal conductivity. At lower 
temperatures, when r(aj = 0) oc Qi^iT /Qy^Y is less than 
T{uj = A) oc eD(A/eD)^ with A ~ V'^/Ey ~ 10"^ ^g. 
ing the small anticrossing induced by the spin-orbit interac- 
tion, the Fermi sea term becomes dominant and the above ap- 
proximation breaks down. As for the longitudinal (Jxx, the fre- 
quency region a; ^ T is relevant and hence 7/2 = r(a; = 0) 
is also justified. Thus, the scaling relation confirmed above is 
expected to hold for T > A even when the structure of the 
self-energy is taken into account. On the other hand, inelas- 
tic scattering by magnons is qualitatively different from that 
by phonons since the electron-magnon interaction involves a 
spin flip. It was shown that the intrinsic contribution is can- 
celed by the side jump contribution and that the skew scatter- 
ing is forbidden in the quasielastic region, i.e., near the Curie 
temperature.^^ 

On the other hand, the intrinsic mechanism is robust against 
inelastic scattering, and its temperature dependence is deter- 
mined only by the Fermi distribution function. Figure 5(a) 
shows that h{T) depends on temperature only under the reso- 
nance condition /i ~ while it is almost independent in the 
case of off resonance. Actually, in the massive Dirac Hamil- 
tonian written as 

Hk = k^a'' - k^'ay + ma\ (10) 

where cr^, cr^, and cr^ are the Pauli matrices, the temperature 
dependence of the Hall conductivity is strong and monotonic 
when the chemical potential is in the gap, but becomes weak 
and nonmonotonic away from the gap, as shown in Fig. 5(b). 
The temperature dependence of the Hall conductivity in fer- 
romagnetic Ni films observed in ref. 19 can be explained by 
the intrinsic mechanism. 

Finally, let us discuss the relevance of our theory to ex- 
periments. In ref. 18, the b term in eq. (3) is responsible for 
the intrinsic mechanism because it is independent of the film 
thickness d. Although we have no parameter corresponding to 
d, h{T) at Xinip = 0.98 is well fitted by ~(Jxy at Ximp = 1 on 
the momentum-space basis in Fig. 5(a), which validates the 
identification of h{T) as the intrinsic mechanism. Therefore, 
our calculations are consistent with the scaling plot in ref. 1 8 
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Fig. 5. (Color online) (a) Temperature dependence of the intrinsic contri- 



bution - 



: 6(T) at /i ~ £2 = 1 foraiin 



: 0.98. Each line is obtained 



by direct calculations of — on the momentum-space basis for Ximp = 1 
and 7 = 0.01. (b) Temperature dependence of the Hall conductivity for the 
massive Dirac Hamiltonian eq. (10). 
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